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ABSTRACT 

An internal shock model is proposed to interpret the radio to infrared (IR) emission of 
the compact jets observed in the hard spectral state of X-ray binaries. Assuming that the 
specific bulk Lorentz factor of the jet at its base varies with a flicker noise power spectrum (i.e. 
P{f) oc 1//), I estimate the energy dissipation profile along the jet and the resulting partially 
self-absorbed synchrotron emission. For this type of velocity fluctuations, and a conical jet 
geometry, the shock dissipation at large distance from the black hole balances exactly the 
adiabatic losses. This leads to a flat radio to IR spectral energy distribution similar to that 
observed in compact jets. 

Key words: accretion, accretion discs - black hole physics - shock waves - relativistic pro- 
cesses - radiation mechanisms: non-thermal - X-rays: binaries 



1 INTRODUCTION 

Internal shocks caused by fluctuations of the outflow velocity are 
, widely believed to power the multi-wavelength emission of jetted 

■ sources such as 7-ray bursts (Rees & Meszaros 1994; Daigne & 
Moscovitch 1998), active galactic nuclei (Rees 1978; Spada et al. 

. 2001), or microquasars (Kaiser, Sunyaev Sl Spruit 2000; lamil et 

■ al. 2010). Internal shocks models usually assume that the jet can be 
discretised into homogeneous ejectas. Those ejectas are injected at 

■ the base of the jet with variable velocities and then propagate along 
" the jet. At some point, the fastest fluctuations start catching up and 
. merging with slower ones. This leads to shocks in which a fraction 

■ of the bulk kinetic velocity of the shells is converted into internal 
energy. Part of the dissipated energy goes into particles accelera- 

„ tion, leading to synchrotron and also, possibly, inverse Compton 
emission. 

Here, we will focus on the applicability of this model to the 
steady compact jets observed in the hard X-ray spectral state of 
black hole and neutron star binaries. These sources have an approx- 
imatively flat Spectral Energy Distribution (SED) extending from 
the radio to the mid-IR (e.g. Fender et al. 2000; Corbel & Fender 
2002; Chaty et al. 2003). These flat spectra are usually ascribed 
to self-absorbed synchrotron emission from conical compact jets 
(Blandford & Konigl 1979) under the assumption of continuous 
energy replenishment of the adiabatic losses. The compensation of 
these energy losses is crucial for maintaining this specific spectral 
shape (Kaiser 2006), and this can possibly be achieved through in- 
ternal shocks. Recently, however, lamil et al. (2010) developed an 
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internal shock model for the emission of jets in X-ray binaries, and 
concluded that energy dissipation only through internal shocks is 
not enough to produce a flat SED. 

In fact this result is not very surprising. Let us consider a jet 
with uniform time-averaged bulk Lorentz factor 7 = 1/-^! — p'^ 
and a section of the jet consisting in a cylindric shell of radius R and 
height H <^ R. This shell propagates downstream, so that at a time 
t, measured in the co-moving frame, its position along the jet is 2: = 
7/3cf. The function R{i) defines the jet geometry. For simplicity, 
in this paper we will consider only conical jets: i.e. R — Rb + 
z tan (j), where (f> is the half-opening angle of the jet, and Rt is the 
radius of the jet at its base. The shell has a specific internal energy e 
(measured in the co-moving frame). The energy losses are usually 
dominated by adiabatic expansion i.e. the pressure work exerted by 
the shell against the external medium as it travels downstream and 
expands: 

dW = PdV = ("ia~l)me-^ ~ — — , (1) 

where P is the total pressure in the shell (including magnetic pres- 
sure), V is its co-moving volume, m its mass and 7a is the effective 
adiabatic index of the flow. For simplicity, it is assumed here that all 
the components of the flow (thermal particles, relativistic particles, 
magnetic field...) have the same adiabatic index corresponding to 
that of a relativistic gas (i.e. 7^ = 4/3). The last term of equation 
[T]is obtained neglecting the effects of possible longitudinal expan- 
sion. At large z, R ~ z tan (f), and, in the comoving frame, the shell 
loses its specific energy at a rate: 

/de\ 2(7„-l)g dfl^ 2g 

\di)ad R di 3i' 
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If, as required to reproduce the flat SED of compact jets, the adia- 
batic losses are fully replenished by some dissipative process, then 
e must be a constant. Equations|2]then implies that the replenishing 
mechanism must have a dissipation profile ex over a long sec- 
tion of the jet. Usual internal shock models predict instead a dissi- 
pation profile scaling like i'^^"^ (Beloborodov 2000) and therefore 
cannot compensate for the losses over a long section of the jet. 

Nevertheless, most studies of the internal shock model so far, 
including that of Beloborodov (2000) and Jamil et al. (2010), have 
implicitly assumed that the Fourier Power Spectral Density (PSD) 
of the velocity fluctuations injected at the base of the jet is flat (i.e. 
white noise). In fact, the energy dissipation profile of the internal 
shocks is very sensitive to the shape of the PSD of the velocity fluc- 
tuations. Indeed, let us consider a fluctuation of the jet velocity of 
amplitude /S.v occuring on a time scale At. This leads to the for- 
mation of a shock at a downstream distance Zs oc At/Au. In this 
shock the fraction of the kinetic energy converted into internal en- 
ergy will be larger for larger Ai;. From these simple considerations 
we see that the distribution of the velocity fluctuation amplitudes 
over their time scales (i.e. the PSD) is going to determine where 
and in which amount the energy of the internal shocks is deposited. 
In this paper, instead of a white noise PSD, we will assume that the 
PSD of the injected fluctuations decreases with Fourier frequency 
like -P(/) oc 1// (i.e. flicker noise). Such noise occurs in many 
physical, biological and economic systems (see e.g. Press 1978). 
In astronomy, it is observed in the solar activity (e.g. Ryabov et al. 
1997) and, most notably, in the X-ray variability of X-ray binaries 
(see e.g. Gilfanov 2010). In Section|2] we show that shocks driven 
by this type of fluctuations can compensate exactly for the adia- 
batic losses in a conical geometry. In Section[3] we estimate the jet 
volume filling factor. Then, in Section [4] we derive estimates for 
the predicted SED of the jet. Finally, in Section|5] these results are 
discussed in the context of Black Hole Binaries (BHBs). 



2 INTERNAL SHOCK DISSIPATION PROFILE 

We assume that small, time dependent fluctuations of the Lorentz 
factor are continuously injected at the base of the jet. The fluctua- 
tions depends on time t and have a Fourrier power density spectrum 
S{f) — Sof^^ for frequencies ranging from /i to /o. The initial 
variance of the injected fluctuations is therefore: 

7LsO= r"2S(/)4f = 2Soln(/o//i). (3) 
Jfi 

We model the jet as a set of discrete homogeneous shells that 
are ejected with such variable velocities at time intervals At — 
(2/o)~^. For simplicity, the mass of the ejectas is assumed to be a 
constant mo = 0.5PjAt/ (7 — l)c'^, where Pj is the total kinetic 
power of the two-sided jet. At injection, the center of momentum 
of two neighbouring shells are separated by a distance Ao — /3cAt. 
This length also corresponds to the smallest scale of the velocity 
fluctuations. As the ejectas propagate downstream, the fastest shells 
catch up and merge with slower ones. During this process of hier- 
archical merging the mass of the ejectas and their separation will 
increase X{t)/Xo = m(t)/mo = K{t). This growth in length 
scale implies a damping of the fluctuations of frequencies higher 
than fo/ K. As a consequence, the variance of the fluctuations de- 
creases: 



Jo 



2 

Trms 



(K) = 



25(/)d/ = 7?msO-25olnA'. 



(4) 



The problem appears simpler when viewed in the frame mov- 
ing with a Lorentz factor 7. Throughout this paper, the tilded quan- 
tities are measured in this co-moving frame. In the limit of small 
scales fluctuations, this frame coincides with the frame of the cen- 
tre of momentum of the shells. In this frame, the average velocity 
of the shells is and its rms amplitude is, to first order in 7rms/7, 



7rms(-R")c 



(5) 



In the limit of low amplitude fluctuations, the velocities in the mov- 
ing frame are non-relativistic. All of the kinetic energy of the shells 
is available for conversion into internal energy. We define the spe- 
cific free energy of the system (i.e. available for dissipation) as 
u{K) — v^^s{K)/2. Combining equations |4] and [5] we see that 
the specific free energy decreases as the scale of the fluctuations 
increases: 



du 

Ik 



(6) 



The fluctuations of scale K collide and merge after travelling dur- 
ing a time: 



i{K) = y\/Av, 



(7) 



where j/ is a factor of the order of unity accounting for the effects 
of the longitudinal extension of the ejecta. It will be estimated in 
section[3] Av is the average (absolute) velocity difference between 
the merging shells. In the case of a linear variability process it can 
be estimated as: 



A ~2 

Aw 



with 



Sqc' 
272/32 



(8) 



fa/K 



8 1 [1 - cos {KTvf/fo)] r^df. 

hi 



G{Kf 

In the limit K <C fo/ fi, the function G{K) is nearly constant: 



(9) 



G{K) ~M = 



-8E 



(2n)!(2n) 



3.6 



(10) 



Then, combining equations |7l[8l[To] we find that the growth of the 
fluctuations is linear: K — t/to, with a characteristic time-scale: 



(11) 



'fi 



to = %/2y7^/3Vo"'-M"'So ^ 

Then combing this with equation|6l we conclude that at times to ^ 
t tofo/fi, the free energy is dissipated at a rate: 

^ = -^.. (12) 

dt 72^2^ 

This scaling of the dissipation profile balances the adiabatic losses 
(equation|2j, for: 

{-Ya - 1) 272/32 

From equation [T3] we can estimate the sound speed in the jet as 
Vs ~ \/ e(7a — 1). Note that, since e also contains the magnetic 
energy, ii^ represents in fact an upper limit on the true sound veloc- 
ity. Using equation [8] and [To] we find the average Mach number of 
the colliding ejectas is at least Avjvs = — 3.6. Such super- 
sonic collisions, ensure that strong shock waves will be generated. 
The region of the jet with significant shock dissipation starts at a 
distance zo — 7/3cto. Therefore, at 2; < 2O, e ~ 0. At the other 
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end of the jet, at distances z > Zf = zofi/ fo, all the energy of 
the internal shocks has been dissipated, and the flow simply cools 
down through adiabatic losses: e ~ ea{z/zf)^^^^ .In the range zo- 
Zf, the specific energy of the flow is approximately uniform and is 
given by equation [T3l 



3 LONGITUDINAL EXTENSION OF THE EJECTAS 

The volume filling factor of the colliding shells, /v, is defined so 
that the average length of a shell of separation A is /vA. If /v <^ 1, 
the jet is constituted of thin colliding 'pancakes', while for /v = 1, 
the jet is a continuous flow. In principle, /v depends on the prop- 
erties of the shells at injection. But /v is also expected to depend 
on z (or t). Indeed, as the shells propagate and gain internal energy 
through shocks, they also expand under the effect of their own pres- 
sure. Here, as it is customary, we assume that this expansion occurs 
at the sound speed in the co-moving frame Vs- For simplicity, we 
will consider that the shells expand in near vacuum, and therefore 
the adiabatic losses due to longitudinal expansion can be neglected 
(free expansion). Any pressure work done during this expansion 
can only be used to accelerate or compress neighbouring shells so 
that overall, this does not represent a loss of energy for the flow. 
We thus do not expect that taking these effects into account would 
change the results dramatically. 

The shells can also undergo significant compression during 
shocks. When two shells are in the process of merging the com- 
pression velocity is simply the difference of velocity of the two 
merging ejectas. During a collision however, the compression lasts 
only for the time necessary for the shock to cross the shell, that is 
~ xfvX/Av, where x~^Av is the speed of the shock in the frame 
of the unshocked ejecta. Applying the jump conditions of Bland- 
ford & McKee (1976) and to first order in Av, their equation (5) 
gives: 



Vl/4 + 7„/2 + (7./2)2 = 7/6. 



(14) 



Since the coUisional time-scale is ~ y\/Av, the shells are com- 
pressed during a fraction ~ xf^/y of the time, and the time- 
averaged compression velocity is reduced to -Ccomp — xf^Av/y 
The parameter y in equationQcan be approximated as: 

-/v 



y = 



(15) 



1 + 

The numerator on the right hand side of equation [Tsl accounts for 
the reduced distance between shells due to their longitudinal ex- 
tension. The denominator corrects for the increased collision rate 
caused by the expansion velocity of the shells. 

As shown in Section |2l the choice of a flicker noise, implies 
that, at large z, the expansion velocity is a constant and indepen- 
dent of /v. The compression velocity instead increases monotoni- 
cally with /v. As a consequence the volume filling factor evolves 
toward an equilibrium in which, on average, compression balances 
longitudinal expansion. This occurs for 

/v ~ [1 + k(1 + >1/2)]"' -0.3, (16) 

independently of the parameters. In fact, the evolution of /v is more 
complex, but in practice, assuming it is a constant turns out to be a 
good approximation. We will therefore set /v = 0.3 and y — 0.5. 



4 SPECTRAL ENERGY DISTRIBUTION 

We now use a simple emission model to estimate the properties 



of the jet SED as a function of the amplitude and bandwidth of 
the Lorentz factor fluctuations. A fraction of the internal energy 
is in the form of relativistic electrons. For simplicity we assume 
that the distribution of the electrons on their Lorentz factors 7e, is 
a power law: n(7e) = no7~^ with in the range 7min- 7max. 
The remaining fraction is in the form of magnetic field energy. The 
particle energy density is: 



nomeC 



7, 



(17) 



where ^ is the equipartition factor that we assume constant. The 
magnetic field energy density is: 

Stt 



1+r 

where the mass density along the jet is given by: 



(18) 



(19) 



2/v(7- I)c27/3c7r' 

As a consequence, in this regime, the magnetic field scales as: B = 
BoRo I R, where _Ro = zorj tan <j), and: 



Bo = 



2foSoM /P,/(7/3)-9(7- 



2/??tan(/i V (7a - l)(l + 0/vC^' 



(20) 



where rj = 1 + Rb / {zo taxi tf)) ~ 1. This scaling holds approxi- 
mately forr > R/Ro > 1, where: r = /o//i + l-»7"^ ^ fo/fi- 
At larger distances, i.e. R > rRo, the magnetic field decays faster 
due to the absence of shock dissipation to balance the adiabatic 
losses. 

The relativistic electrons emit synchrotron radiation. Their 
emissivity at a co-moving photon frequency D is (Rybicki & Light- 
man 1979): 



P+5 p~ 1 

j,^Kj^B—i>- — .^ 



where 



p/ 3p+19 \ ^3p-l 



12 



12 



(21) 



(22) 



and r is the usual gamma function. This emission can be self- 
absorbed. The absorption coefficient is : 



(23) 



/ 3e / 3p + 2 \ / 3p + 22 \ 

\2TrmecJ V 12 / V 12 y 



The specific intensity at the surface of the jet is: 



.(24) 



(25) 



In the case of a conical jet in a steady state, making an angle 6 with 
the line of sight, the flux at an observed frequency u — Sv, is: 



tan < 



Rh{v)dR, 



(26) 



where 5 — [7(1 — /3cos6')]~^ is the Doppler factor and D the 
distance to the source. In order to obtain analytical estimates of the 
observed SED we now assume that the magnetic field has the R~^ 
dependence given by equation |20] for R between _Ro and Rf — 
rRo, and B = elsewhere. The SED can then be approximated 

as follows: At the highest frequencies the emission is optically thin 

_ p-i 

and the spectrum is a power law oc u 2 , while at the lowest 
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10'' 10" 10'° 10" 10'' 10'' 10'-* 10'=^ 

V [Hz] 

Figure 1. Typical SED obtained for a source at a distance of 2 kpc and an inclination of 40° . The dashed curve shows the result of a Monte-Carlo simulation, 
the full lines shows the analytical approximation. The jet lorentz factor is 7 = 2 for a jet kinetic power Pj = lO^^L^ fluctuating with a fractional rms 
ampHtude of 30 percent with a flicker noise PSD at Fourrier frequencies ranging from /i = 1.6 X 10~^ Hz to /o = 50 Hz. The jet radius at the base is 
Rb = 10/?G 2nd its half opening angle is = 1°. The equipartition ratio is ^ = 1, the electrons in the jet have a power-law energy distribution with slope 
p = 2.3 between 7niin = 1 and 7inax = 10®. The emission from the counter jet is not shown. 



frequencies the emission is optically thick Fu cc . Forp > —3, 
there is a range of intermediate frequencies, comprised between Vs 
and Vt, in which the emission at a given frequency is a mixture of 
optically thick and thin emission from different regions of the jets. 
In this regime, the flux is independent of the photon frequency (see 
e.g. Kaiser 2006) : 



(27) 



where a = —5/{p + 4). The observed turnover frequency, vt, is at 
the intersection of the optically thin and flat asymptotic branches: 



I't 



5 (K^iRoBt'^ 



2 



(a+ l)r(a) 



(28) 



The observed frequency of transition to the low frequency optically 
thick regime is given by: 



S K^iRoB'^ 



2 

\ p + 4 


r2 — 1 




aV{a) 



(29) 



These analytical estimates have been tested against Monte-Carlo 
simulations of the internal shock model with a code similar to that 
of Jamil et al. (2010). Fig. 1 compares the average SED given by the 
analytical formulae to that obtained from a simulation with param- 
eters corresponding to Cyg X-1 (see caption and Section [5}. The 
main differences occur around and below Ps when the self-similar 
approximation that we used in the analytical formulae breaks down. 
Also, unlike the simulation, the analytical model ignores the emis- 
sion of from the jet at distance larger than z / . 

It is worth noting, however, that in order to obtain analytical 
estimates of the SED, we had to make several simplifications. First, 
equation [25] neglects the angular dependence of the optical depth, 
this is a good approximation only for systems with a large inclina- 
tion. Then, we also assumed that the electron energy distribution 
is a pure power-law and therefore neglected the effects of a ther- 
mal component in the electron distribution. We also neglected the 
effects of adiabatic cooling on the shape of the electron distribu- 
tion. Simple estimates indicate that both effects will affect the SED 
only marginally. On the other hand the effects of radiative cooling 
can be very strong in the region close to the base of the emitting 
region (see e.g. Chaty et al. 2011). These radiation losses, that we 
have neglected, affect predominantly the most energetics electrons. 



emitting in the optically thin regime, while the flat optically thick 
part of the SED is essentially unchanged (Zdziarski et al. 2012). 



5 APPLICATION TO BLACK HOLE BINARIES 

We can now use the formalism developed in the previous section to 
estimate the the jet parameters and emission properties in BHBs. In 
the following we assume a black hole of mass M = lOmiM© and 
the radius of the jet at 2: = is: 



Rb 



IQrxRa = 1.5 X lO^riTTii cm. 



(30) 



If the jet is launched by the accretion disc, the highest and lowest 
frequencies of variability /o and /i must relate to time scales of the 
accretion flow. The picture could be similar to that of the model by 
Lyubarski (1995) in which density fluctuations generated at a range 
of disc radii propagate inward and lead to flicker noise fluctuations 
of the emission of the inner part of the accretion flow. Similar fluc- 
tuations may also lead to a variable velocity of the ejected shells. 
The detail of how the disc fluctuations are transferred to the jet is 
out of the scope of this paper. An interesting possibility is that the 
rapid episodic ejections could be produced through MHD mecha- 
nisms similar to that leading to coronal mass ejections in the sun 
(Yuan et al. 2009). Let g be the ratio of the ejection time scale 
At — /(7^/2 to the Keplerian time-scale of the accretion disc at 
radius Rb- In our numerical estimates we will assume that those 
time-scales are comparable (i.e. g ~ 1): 



/o ~ 50 ' m 



3/2 -1 -1 „ 
' " Hz. 



(31) 



Similarly we associate the lowest frequency of the fluctuations to 
2g times the dynamical time-scale at a large, possibly outer, radius 
of the disc Rd'. 



Rd = lO^rs-Rc = 1-5 X 10^ rgmi cm, 
/i ~ 5 X 10"^ r'^^'^m^^g'^ Hz. 



(32) 
(33) 



We will assume that the fractional rms amplitude of fluctuations of 
the kinetic power is a = 7rmsO /(7 ~ 1) ~ 0.3, i.e. comparable to 
that observed in the X-ray variability of black hole binaries in the 
hard state. The average travel time of a shell from the ejection point 
to the first shocks is then: 



7 + 1 gr?''^miC^/^ 
to ^90 ^^-t ms, 

3 (70.3 



(34) 
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where (T0.3 — (t/0.3, and C, = In (.Rd/i?!,)/ In 10*. This time- 
scale might be related to the delays of comparable amplitude ob- 
served between the X-ray emission of the accretion flow and the 
optical and IR emission from the jet in GX339-4 (Casella et al. 
2010; Gandhi et al. 2010) or XTE Jl 1 18-H480 (Kanbach et al. 2001 ; 
Hynes et al. 2003; Malzac et al. 2003). It corresponds to a lowest 
distance at which significant emission takes place of: 



2:0 ~ 5 X 10' 



9 (7 + 1)7/? grTm.e'^ 



27 



O"0.3 



(35) 



This is comparable to the distance of ~ 10 -Rg inferred in Cyg X-1 
both by the modelling the of SED (Zdziarski et al. 2012), and that 
of the orbital modulation of the radio emission (Zdziarski 2012). 
The size of the dissipation region is then 



z/ ~ 5 X 10' 



(7 + 1)7/3 grTrn^e'^ 



27 



O"0.3 



(36) 



which is roughly comparable to the extension of the radio jet of 
Cyg X-1 in the VLBA images of Stirling et al. (2001). These au- 
thors also constrain the jet opening angle of Cyg X-1 to be small 
</> = '/'lyfo < 2°. From this, we can estimate the radius of the jet 
at the base of the emitting region: 



Ro 



-8x10'^ (7+1)7/3 <l)nigrl''^m,iC}''^ 
V27 



O"0.3 



(37) 



Ha and [Olll] measurements of the optical nebula surrounding 
Cyg X- 1 (Gallo et al. 2006 ; Russell et al. 2007) indicate that the 
jet power is of the oder of a few percent of the Eddington luminos- 
ity Le- This implies a magnetic field at the base of the emitting 
region: 



Bo~2x 10 



where Pv, 



„2 
'^0.3 

01 



162^ (7 + 1)" 



% 



1 + i 7^/3^ 



r?mi 



1/2 



(38) 



— ^ oii-E ■ '^'^^^^ estimates for Ro and So are compa- 
rable to the values inferred in Cyg X-1 (Zdziarski et al. 2012), GX 
339-4 (Gandhi et al. 201 1) and XTE J1550-564 (Chaty et al. 2011). 
Fixing the parameters of the electron distribution to standard values 
{p = 2.3, 7niin ~ 1, 7max ~ 10^) , we Can now evaluate the flux 
of the flat section of the spectrum: 



2^/27 P%mi(7^.3 



1 + CC(7 + 1)7/3 



2p+13 
2p+fi 



mjy, 



(39) 



where Dkpc is the distance to the source expressed in kpc. This 
estimate is consistent with the flux observed in bright hard states 
(e.g. at ~ 2 kpc, Cyg X-1 has an average radio flux of of 15 mJ at 
15 GHz). Moreover we see that for p ~ 2 the flux depends on the 
jet power approximately as F,^o oc P} *, as inferred from the ob- 
servations of the radio-X-ray correlation (Gallo, Fender & Pooley 
2003). The turnover frequency is expected in the mid-infrared, in 
agreement with the observations of Cyg X-1 (Rahoui et al. 2011) 
or GX 339-4 (Gandhi et al. 201 1): 



ft ~ 1.8 X 10" mr, 



where 



m = 



2p + 10 
P + 4 



p + 5 



Hz, 



T + I 



(40) 



P+6 

2p+S 



3^ 



3p + 14 
2p+8 



(7+1)7/3 



.(41) 



Finally the model predicts that the emission should decrease signif- 
icantly at frequencies below: 



(42) 

and this is something that can be investigated in the near future with 
LOFAR. Note that in Fig. 1, i/s occurs at much higher frequencies 
(a few GHz). Indeed, because simulating very long time-scale fluc- 
tuations is time consuming, fi was set to 1.6 x 10^'^ Hz in Fig. 1, 
rather than ~ 5 x 10^^ Hz as given by equation[33l 

We conclude that internal shocks driven by flicker noise fluc- 
tuations can produce not only the flat SED of compact jets in BHBs, 
but also other properties such as the flux amplitude or the location 
of the break frequency. The model also predicts strong multiwave- 
length variability that will be the focus of future works. 
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